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1 .  Introduction 

Currently,  about  25  percent  of  the  energy  consumed  in  the  United 
States  is  in  the  form  of  electric  energy  generated  by  power  plants.  More 
than  80  percent  of  this  electricity  production  is  due  to  steam  power  plants 
[15] .  This  reveals  the  important  role  of  steam  turbine  power  plants  as  a 
major  contributor  to  energy  production  and  suggests  the  efficient  design 
and  operations  of  these  plants  as  a  vitally  important  objective  in  this 
era  of  energy  shortage. 

The  efficient  design  of  steam  power  plants  has  been  a  source  of 
contention  since  the  late  19th  century.  During  the  past  80  years,  the 
efficiency  of  the  steam  power  plants  has  been  increased  from  about  4  per¬ 
cent  to  about  33  percent.  That  is,  only  about  one-eight  as  much  fossil 
fuel  is  now  required  per  unit  of  electric  energy  production,  compared  to 
that  of  1900  [12J.  This  advance  has  been  due  to  a  variety  of  reasons, 
such  as  mechanical  design  improvements,  advances  of  thermodynamics,  and 
metallurgical  developments  that  have  made  high  temperature  boilers  and 
turbines  possible. 
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There  liave  been  numerous  studies  addressed  to  the  optimal  design 
and  economic  operation  of  the  various  components  of  steam  turbine  power 
plants.  A  good  survey  of  these  studies  is  given  by  Wiebking  [15].  The 
major  components  consist  of  the  boiler,  steam  turbine,  condenser,  pump 
and  generator.  Figure  1  depicts  a  simpiified  layout  of  these  components 
in  a  steam  turbine  power  plant. 


Figure  1 .  A  simplif ied  layout  of  a  steam 
turbine  power  plant 

It  is  welL  known  that  the  optimal  design  of  the  subsystems  and 
individual  components  of  a  given  system  do  not  necessarily  lead  to  an 
overall  optimal  design  of  the  given  system.  In  power  plants,  however, 
Wilson  and  Malouf  [17]  and  Wilson  [161  have  shown  that  the  overall 
optimal  sizing  of  the  condensing  system  and  turbine  exhaust  annulus  can 
be  closely  approximated  by  optimizing  these  systems  independently  of 
the  remaining  components.  Based  on  this  finding  Wiebking  [15]  has 
formulated  and  solved  numerous  geometric  programming  models  with 
different  design  configurations,  for  optimal  sizing  of  the  turbine 
exhaust  and  condensing  systems  of  steam  turbine  power  plants. 

In  che  present  paper,  we  analyze  one  of  these  GP  models,  studied 
by  Ecker  and  Wiebking  [3],  [4],  for  optimal  value  function  sensitivity 
and  parametric  bounds,  with  respect  to  many  of  the  model  parameters. 

The  sensitivity  analysis  methodology  used  here  is  based  on  the  work  of 
Fiacro  [6],  and  the  optimal  value  bound  approach  was  proposed  by  Fiacco 
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and  implemented  by  Fiacco  and  Chaeml  [10].  Ghaemi  has  successfully 
demonstrated  the  applicability  of  this  bound  calculation  procedure  to 
several  nontrivial  problems  [11]. 

It  is  not  the  purpose  of  this  study  to  elaborate  on  and  justify 
the  choice  of  constraints  or  objective  function,  or  their  derivation, 
involved  in  the  model  under  consideration.  Readers  interested  in  these 
aspects  may  refer  to  [15].  Here, after  presenting  the  mathematical  formula¬ 
tion  of  the  model  and  identifying  the  variables  and  parameters,  we  report 
on  the  solution  of  a  convex  equivalent  of  the  model  and  analyze  it  for 
sensitivities  and  optimal  value  bounds,  using  the  latest  version  of 
our  penalty  function-sensitivity  analysis  computer  program  SENSUMT  [10]. 

The  paper  is  organized  as  follows.  Section  2  briefly  reviews 
the  basic  idea  of  the  optimal  value  bound  calculation,  its  implementa¬ 
tion  on  the  computer,  and  its  application  to  GP  models.  Section  3 
lists  the  model  parameters  and  variables  and  problem  formulation.  The 
usual  exponential  transformation  is  applied  to  convert  the  GP  to  a 
standard  convex  nonlinear  programming  (NLP)  problem.  Sections  4  and  5 
present  the  solution,  sensitivity,  and  bound  analysis  results. 

2.  Parametric  Optimal  Value  Bounds 

In  this  section  we  briefly  outline  the  procedure  for  calculating 
piecewise  linear  parametric  upper  and  lower  bounds  on  the  optimal  value 
function  of  the  convex  right-hand  side  perturbation  problem,  which  is 
known  to  have  a  convex  optimal  value  function.  Making  additional  use 
of  some  recent  developments  by  Dembo  [2],  we  will  demonstrate  the 
applicability  of  this  technique  to  geometric  programs. 

Parametric  Bounds  on  the  Optimal  Value  Function  of 
Convex  Right-Hand  Side  Problems  CR(e) 

Consider  the  following  parametric  programming  problem  R(e) : 

Minimize  f(x) 

subject  to  g^x)  £  ei  •  i=l,m  R(e) 

h. j  ( x)  -  ,  j  =rrt-l ,  p  , 
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where  the  decision  variable  xeE  ,  f  ,  g^,  and  :  E  -*■  E  are  C~ 

and  the  parameter  vector  e  =  (e^  ,  ....  ep)  is  in  E^  .  If  £(x)  and 

-g^(x)  »  i*l»m  ,  are  convex  and  h^(x)  ,  j=nrfl,p  are  affine,  then  the 

problem  R(c)  is  convex  and  will  be  designated  as  CR(e)  .  It  is  well 
known  that  f*(e)  ,  the  optimal  value  function  of  the  problem  CR(e) 
is  convex  over  any  convex  subset  of  the  parameter  space  on  which  a 
solution  is  defined. 

The  convexity  of  the  optimal  value  function  f*(e)  of  the  problem 
CR(e)  allows  for  a  simple  calculation  of  parametric  linear  upper  and 
lower  bounds  of  this  function  when  any  of  the  problem  parameters  is  per¬ 
turbed.  In  our  application,  this  will  require  the  solution  and  correspon¬ 
ding  optimal  value  function  sensitivity  information  for  both  perturbed  and 
unperturbed  problems.  The  idea  follows  immediately  from  two  well  known 
properties  of  convex  functions: 

(I)  Any  line  connecting  two  points  on  the  graph  of  a 
convex  function  does  not  underestimate  that 
function  between  the  points. 

(ii)  Any  tangent  line  to  the  graph  of  a  convex  function 
does  not  overestimate  that  function. 

These  two  properties  lend  themselves  in  a  natural  way  to  the  calculation 
of  parametric  bounds  on  the  optimal  value  function  of  the  problem  CR(c) 
under  large  perturbations  of  any  of  the  problem  parameters,  say  . 

The  idea  of  using  properties  (i)  and  (ii)  for  calculating  bounds 
on  the  optimal  value  function  of  the  problem  with  a  convex  or  concave 
optimal  value  function  and  for  estimating  bounds  on  nonconvex  problems 
by  way  of  estimating  problems  with  convex  or  concave  optimal  value  function 
is  due  to  Fiacco  [7].  Fiacco  also  notes  that  techniques  for  generating 
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underestimating  problems  that  are  jointly  convex  in  the  decision 
variable  and  the  parameter  (hence,  that  have  a  convex  optimal  value 
function  113])  are  already  available  for  problems  that  are  jointly 
separable  [5]  or  jointly  factorable  [14],  thus  making  it  possible  to 
calculate  lower  bounds  using  known  procedures.  The  first  computational 
implementation  of  this  approach  was  developed  and  reported  by  Ghaemi  [11], 
who  also  derived  the  formulas  for  calculating  a  jointly  convex  over¬ 
estimating  problem  of  a  factorable  program,  thus  making  it  possible  to 
numerically  implement  the  proposed  scheme  for  calculating  upper  bounds 
when  the  problem  functions  are  factorable  (as  they  inevitably  are,  in 
practice) . 

Suppose  f*(e)  is  differentiable.  (Well-known  conditions  guaran- 
eeing  this  were  given  by  Fiacco  [6].  ) 

Basis  Cor  Calculation  of  Bounds  on  f*(e)  of  the 
Problem  CR(e) 

For. simplicity,  suppose  that  some  e  is  perturbed  from 
(  .  to  t..  ,  while  the  remaining  parameters  are  fixed  at  their 

base  vaiues. 


Step  1.  Solve  the  unperturbed  problem  and  obtain  f * (t •  ) 


and  df*(ei)/dei  at  e±  =  .  Under  general 

conditions  (e.g.,  see  Armacost  and  Fiacco  [1]), 


df*^)  ui(ei>  * 

dGi  ~wi^ep  •  • 
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where  u  and  w  are  optimal  Lagrange  multipliers  associated  with 
inequality  and  equality  constraints  of  CR(e)  ,  respectively. 

Step  2.  Resolve  the  problem  and  obtain  f*(e^)  and 
df*(ei)/dti  at  <  1  =  l12  . 

Step  3.  Derive  the  equation  of  the  line  passing  through  the 

points  (cn  ,  fMe^))  and  (c12  ,  f*(ei2))  -  This 

line  provides  a  parametric  upper  bound  f*(e^)  for 

f*(e.)  as  a  function  of  r.  r  [e.„  ,  e  . 
i  x  il  i2 


Step  4. 


Derive  the  equation  of  the  tangent  lines  to  f*(e.) 
at  the  above  two  points  with  the  slopes 


df*(e.) 

de^ 


and 


P  ~c 

i  il 


df*(cA) 

dc~! 

l 


ei=Ci2 


calculated  in  Steps  1  and  2,  respectively.  The 
maximum  of  these  two  lines  provides  a  piecevise- 
linear  parametric  lower  bound  f*(Lj  for  f*(ej 

as  a  function  of  e.  in  [c.,  ,  e.„l  . 

l  il  i2 


The  lines  obtained  in  Steps  3  and  4  provide  the  desired  bounds, 
forming  a  triangle  which  encloses  the  optimal  value  function  f*(e^) 

over  the  given  range  of  .  Further,  a  smooth  estimate  of  f*(ei) 

over  the  given  interval  can  also  be  made  by  fitting  any  differentiable 
convex  function  that  passes  through  the  points  given  in  Step  3,  having 
the  corresponding  slopes  at  these  points  obtained  in  Steps  1  and  2. 
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In  our  implementations,  we  calculate  the  highest  quadratic  function  pass¬ 
ing  through  the  given  end  points  and  having  the  required  slope  at  one  of 
the  end  points.  It  can  be  shown  that  this  function  is  entirely  contained 
in  the  indicated  triangle  and  is  convex. 

It  is  clear  that  the  fundamental  property  exploited  in  calculating 
the  above  bounds  is  the  convexity  of  the  optimal  value  function.  Thus, 
the  use  of  this  procedure  is  not  limited  to  problems  of  the  form  CR(e), 
but  can  be  applied  to  any  parametric  problem  that  has  a  convex  optimal 
value.  In  fact,  it  is  obvious  that  this  technique  is  also  applicable 
to  parametric  programs  with  concave  optimal  value  functions.  The  only 
minor  alteration  is  that  Step  3  of  the  algorithm  will  provide  the  lower 
bound,  while  Step  4  will  yield  the  upper  bound. 

The  relevance  of  the  application  of  the  above  technique  to  the 
problem  under  consideration  in  the  present  study  is  that  the  posynomial 
programs,  under  the  usual  exponential  transformation,  becomes  a  convex 
problem.  Since  the  transformation  does  not  affect  the  value  of  the 
problem  functions,  the  above  technique  may  obviously  be  applied  to 
calculate  parametric  bounds  on  the  optimal  value  function  of  GP  problems 

with  respect  to  right-hand  side  perturbations.  The  convexity  of  the 

optimal  value  of  the  posynomial  GP  as  a  function  of  the  right-hand  side 

perturbations  follows  immediately  from  the  exponential  transformation 

and  a  well-known  result.  Furthermore,  it  follows  that  the  optimal 

value  function  of  the  primal  GP  (posynomial)  problem  (formulated  as  a 

minimization  problem)  is  a  monotone  nondecreasing  concave  function  of 

the  coefficients  appearing  in  the  GP  primal  objective  function.  This 

follows  from  the  linearity  of  the  objective  function  as  a  function  of 

these  coefficients  and  the  fact  that  f*(t)  =  min  f(x,e)  is  concave 

K 

if  f(x,e)  is  concave  in  e  and  R  does  not  depend  on  z  .  Concavity 
with  respect  to  each  was  noted  by  Derabo  [2],  who  did  not  make  use 

of  the  general  result  indicated  but  calculated  the  second  derivatives  of 
the  optimal  value  function.  In  Section  4,  using  these  results,  we  are 
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able  co  apply  the  given  technique  to  derive  piecewise  linear  parametric 
bounds  on  the  optimal  value  function  of  the  model  under  consideration  as 
a  function  of  the  model  parameters. 

3.  The  Model 

As  mentioned  earlier,  it  is  not  the  purpose  of  this  study  to 
justify  the  underlying  theory  and  formulation  of  the  resulting  GP  problem 
We  list  the  problem  variables  and  parameters,  followed  by  a  concise  state¬ 
ment  of  the  model  objective  function  and  constraints.  Readers  interested 
in  details  regarding  the  derivation  of  the  model  may  refer  to  [3],  [4] 

and  [15]. 

The  model  under  consideration  requires  the  minimization  of  the 
total  cost  involved  in  the  steam  turbine  power  plant  as  a  function  of  the 
turbine  exhaust  annulus  and  the  condenser  system  design,  subject  to  a 
multitude  of  engineering  and  thermodynamic  constraints. 

Model  variables 

The  variables  involved  in  the  model  and  their  description, 

extracted  from  [3],  are  as  follows: 

A  :  total  turbine  exhaust  annulus  area,  sq  ft 
an 

T  :  saturation  temperature  of  steam, °F 
s 

N  :  number  of  tubes  in  condenser 

D  :  outside  tube  diameter,  ft 

o 

D,  :  inside  tube  diameter,  ft 
L  :  condenser  tube  length,  ft 

q  :  rate  of  heat  transfer,  Btu/hr 

y  :  condenser  flow, lb/hr 

cl 

:  moisture  correction  factor  at  expansion  line 
end  point 


-  8  - 


« 


T-4  37 


Ah  :  exhaust  loss,  Btu/lb 

:  hot  water  temperature,  °F 

Model  parameters 

There  are  numerous  parameters  involved  in  the  derivations  of  the 
models  developed  in  [15]  by  Wiebking.  Here  we  will  only  state  those 
parameters  which  are  explicitly  involved  in  the  formulation  of  the  model 
considered  here: 


a^(i=l,6):  regression  coefficient  involved  in 

derivation  of  emperical  formulas  for 
the  model 


b . (1=1,5  : 
1  and 
1=61, 
62) 


an 


an 


D 


o,  ciin 


M,R 


2,R 


k 

£ 

Pr 


pressure  drop  factor 

regression  coefficients  involved  in 
derivation  of  emperical  formulas  for 
the  model 

exponent  in  annulus  area  cost  relation 

coefficient  in  annulus  area  cost  relation 

cost  of  electricity,  $/(kW-hr) 

specific  heat,Btu/(lb-°F) 

unit  cost  of  condenser  surface,  $/sq  ft 

minimum  allowable  outside  tube  diameter,  ft 

moisture  correction  factor  at  rated  average 
conditions 

2 

gravitational  constant,  ft /sec 

expansion  line,  end  point,  at  rated  average 
conditions,  Btu/lb 

thermal  conductivity,  Btu/f t-hr-°F) 

tube  wall  thickness,  ft 

Prandtl  number 
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P^  :  rate  of  depreciation  per  year 

Q  :  condenser  flow  aL  rated  average  conditions,  lb/hr 
a ,  K 

:  cold  water  temperature,  °F 

t  :  load  factor  -  effective  operating  time,  hours/year 

u  .  :  minimum  tube  velocity,  ft /sec 

min 

V|  :  lower  bound  on  annulus  velocity  of  steam,  ft/sec 

Vy  :  upper  bound  on  annulus  velocity  of  steam,  ft/sec 

W  :  water  flow  rate,  lb/hr 

Hg  :  generator  efficiency 

Hp  :  pump  efficiency 

A  :  latent  heat  of  evaporation,  Btu/lb 

M  :  viscosity,  lb/(ft-hr) 

p  :  density,  ib/cu  ft 

Tmax  :  allowable  last  stage  bucket  loading,  lb/(hr-sq  ft) 

^fmax*  maximuni0aAlowable  temperature  increase  of  cooiing 
water,  r 

f  :  outside  film,  a  subscript 

The  objective  function  gQ 

The  model  minimixes  g^  ,  the  total  annual  cost  involved  in  the 

steam  turbine  power  plant  as  a  function  of  the  exhaust  annulus  and  the 
condenser  system  design,  and  consists  of  two  parts:  fixed  costs  and 
operating  costs.  Denoting  the  fixed  and  operating  costs  of  the  turbine 
by  Cp  <r  and  C0  and  those  of  the  condenser  by  Cp  ^  and  CQ  ^  , 

respectively,  the  total  annual  cost  to  be  minimized  is 
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The  formulas  derived  in  [31  and  [A]  for  the  above  cost  components, 
as  a  function  of  the  design  variables  anil  parameters,  are  as  follows: 


:F,T 

C1 

A  15 

A  an 

an 

$/year 

0,  T 

C2 

(q-q^)  +  fuel  cost 

$/ year 

^.C 

C3 

D  LN 

0 

$/year 

'0,c 

C4 

L/fD^-8  *  N1'8) 

$/year 

where 

c,  =  c  P 
1  an  c 

c 2  =  fpt  Hg/ 3412 . 75 

c  =  tic  P 
3  t  c 

c4  =  (32)  ( .046)  CEBt  W2-8  (p/A)' 2/(gnp  p2  tr1’8) 

and  q  is  the  rate  of  heat  transfer  at  the  rated  average  conditions 
K 

in  Btu/hr. 

Notice  that  the  operation  cost  of  the  turbine  Cq  consists  of 

the  replacement  energy  cost  and  the  fuel  cost.  In  the  formulation  of 
the  mode],  it  is  assumed  that  the  turbine  generator  runs  at  full  load, 
implying  that  the  cost  of  fuel  is  constant.  Also,  in  replacement  cost 
terms,  the  quantity  c^q^  is  constant.  Thus,  the  fuel  cost  and 

are  deleted  from  consideration  in  the  optimation  model. 

The  Constraints 


(i)  constraint  on  condenser  rejected  heat 


T- 


whoro 


s 


■>2.K  +  Vh.R  +  V*  -  “2 


Vm.R  +  Vo  -  b2 


c;  -1 


Ah*  =  Ah  -  a. 


where 


(ii)  i4oisture  correction  factor  constraint 

.  T 

1  s  ,  , 

82  :  C8  ”  C9  ”1  -  1 

1  M  M 


c„  =  .87  a. 

8  4 


c9  =  *87  b4 


(iii)  Exhaust  loss  constraint 

-  2b s 


83  Cl°  A2  Ah’ 
an 


+  c 


11  b 

T  Ah’ 


<  1 


where 


10 


'11 


*  b 


62 


a5  Qa,R 


(iv)  Heat  exchanger  mean  temperature  drop  constraint 

4/3  g  D. 


g,.  :  c 


12  T 


+  c 


13  T 


C14  _  „7/6,.  ,4/3  +c15_.:.-2 


T  N  '  D  L 
s  o 


T  L  N 
s 


-  12  - 


where 


where 


where 


T-4 


12 


13 


'14 


15 


1/(2  cp  W) 


(.725  *) 


-4/3 


4/3 


(2) (3600) 2  kf3p  f2g  X 


-1  .8 

(.023  *  k  Pr  )  *  (^) 


(v)  Condenser  flow  constraint 


55  ' 


16 


+  c.  _  T  <  1 

17  s  - 


C16 

C17 


and 


(vi)  Heat  exchanger  tube  diameter  constraints 


% 


8 


7 


c18 

C20 


19 


1 

D 

o 


<  1 


<  1  and 


21 

1. 


'20 


o,min 
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(vii)  Annulus  fluid  velocity  constraints 


where 


where 


where 


&8  :  c2 1 


<  1 


an 


g9  :  c22 


an 


<  1 


"21 


"22 


ac  Q 
5  va,R 


(viii)  Last  stage  bucket  loading  constraint 


s10  •  c23 


<  1 


an 


c 


2  3 


1 


T 

max 


<i*>  Heat  exchanger  tube  water  velocity  constraint 
811  :  c24  “  J  1 


c 


24 


3600  u  J 
min 

4W 


TT  p 
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(x)  Thermal  pollution  constraint 


*12  =  C25  Th  i  1 


T  +  AT 
c  max 


Denoting 


j  the  model  variables  Afln  ,  Tg  ,  N  ,  Do  ,  Dj  , 


L  ,  q  ,  Qa  ,  fM  .  <'&  %>  Th  b>  xi  through  xu  , 

respectively,  the  mode!  4*  i&i  iiulated  as  follows. 

The  GP  problem 


-A. 8  -1.8 


Minimize  g„  *  0,,  «x  '  +  «2  *7  +  C3  x3  x4  X6  +  C4  x6  X5  X3 

x  1  *1 

subject  to:  Bl  =  c5  xg  x~  +  cg  xg  x2  x?  +  c?  xg  x9  x1Q  x?  < 

g2  =  c8  x9  1  +  c9  x2  x9  <1 

2b5  _2  -1  "b5  -1 

«3  '  C10  X2  X1  X10  +  CU  X1  X2  X10 

_i  -1  4/3  -1  -7/6 

g4  -  cl2  x2  +  cn  x?  x2  +  x7  x2  x3 

.8  -1  -1  -.2  <  1 
x  -1  x  -4/3  +  c15  xy  x5  x2  x6  x3 

4  6 


c,,  x„  +  c^7  x2  <  1 


g5  =  C16  X8 


*6  2  c18  X4_1  +  C19  X5  X4_1  2  1 


g7  =  C20  X4  -  1 


°5  -l 

*8  "  C21  X2  X1 


g9  -  c22  xL  x2 


-  15  - 


V  | 


T-43? 


810 


23  X8  "1 


-1 


1 


811  *24  X3  X52 


<  1 


812 


C25  X11 


< 


1 


4.  Model  Solution 


Ecker  and  Wlebklng  [3],  using  a  data  base  drawn  from  [15] ,  have 
solved  this  model.  In  [3],  they  conclude  from  a  preliminary  analysis 
that  constraints  gg  -  g^2  i  with  the  given  data  base,  are  not  binding 

and  thus  ignore  these  constraints  in  determining  a  solution  to  the  CP 

problem.  Henceforward,  the  problem  with  gQ  -  g  deleted  will  be  the 

o  1 Z 

focus  of  our  attention,  and  will  be  called  the  "reduced  problem." 


The  values  of  coefficients  c^  through  c2Q  in  [3]  are  given 
as  follows: 


2541.0 

C6  = 

.23972 

012293 

C7  = 

1.0 

1 

2.4190 

C8  “ 

.62004 

77171.0 

C9  = 

1.1072  * 

lO"3 

888.76 

C10= 

2.1872  * 

ioife 

C12 
C13 
C14 
c]  5 


=  6.2139  *  10 

=  50.0 

=  4.7394  *  10' 
=  7.3124  *  10 
=  1.2577  *  10‘ 


c11L  =  2.447  *  106 
lo 

Cjl7  =  1.0289  *  10 
clg  =  8.1667  *  10 
cJ9=i.O 
c2Q  =  .083333 


Using  the  above  coefficients,  they  obtain  the  solution  of  the  reduced 
model  by  means  of  a  modified  concave  simplex  algorithm  that  solves  the  dual 
geometric  program.  The  solution,  obtained  after  300  seconds  of  CPU  time 
on  an  IBM  system  360/50  computer,  is  listed  below. 

Dual  objective  function:  32.3728  ■  10^ 
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Pr imal 

objective  function: 

32.  oj5b  •  101' 

Primal 

variables: 

q 

2.5503  •  109 

T  =  85.455 
s 

A 

an 

409.16 

N  =  13,662 

fN  " 

0.71520 

D.  =  0.082461 

l 

= 

2.7010  •  106 

D  =  0.091950 

0 

Ah'  =* 

46.882 

L  =  70.606 

To  solve  the  reduced  GP  problem  in  this  study,  the  convex  equiva¬ 
lent  of  the  (reduced)  model  was  formulated  via  the  usual  exponential 

-t . 

transformation,  x.  =  e  ,  i=l  ,  10  ,  and  the  same  coefficients  c. 

l  *  ’  l 

used  by  Ecker  and  Wiebking  and  reported  in  [3]  were  utilized.  To  prevent 
over-and  under-flow  in  the  process  of  obtaining  the  numerical  solution, 
the  following  constraints  (which  proved  to  be  nonbinding  at  the  optimal 
solution)  were  added  to  the  model: 

-  170  <  t , -4 . 8  tc  -  1.8  t„  <  170 

-  o  5  3  - 

-  3°  <  t  <  30  i=l, 10 


The  resulting  problem  was  solved  by  SENSUMT,  on  The  George  Washington 
University  IBM  3031.  The  solution,  obtained  in  57  seconds  of  CPU  teme,  is 
listed  below. 


Dual  objective  function 
Primal  objective  function 


32.37324  *  10 


32.37325  *  10 


X1 

=  A 

an 

=  408.88 

x6  = 

L  = 

58.46230 

X2 

=  T 

s 

=  85.452 

X7  “ 

q 

=  2.52929  * 

10 

X3 

=  N 

=  18979.5 

x8 

Qa 

=  2.682*8  * 

10 

x4 

=  1) 

o 

=  .07902 

x9  * 

fM 

-  .  71 455 

x5 

“  °i 

=  .07085 

X10  ” 

Ah’ 

=  46.88288 

9 

6 
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The  significantly  greater  accuracy  of  this  solution,  based  on  the 
gap  between  the  values  of  the  dual  and  primal  objective  functions,  compared 
to  that  obtained  in  [3],  is  apparent  from  the  given  data.  The  two  solution 
vectors,  although  considerably  different  in  some  components,  are  of  the 
same  order  of  magnitude.  As  noted,  the  CPU  time  required  by  SENSUMT  to 
obtain  a  relatively  very  accurate  solution  of  the  given  model  was  about 
1/5  of  the  CPU  time  spent  by  the  algorithm  in  (3]  to  obtain  a  rather  crude 
estimate  of  a  solution.  Of  course,  since  the  IBM  3031,  is  faster  than  the 
IBM  360/50  and  since  we  do  not  know  the  initial  point  used  in  [3]  to  solve 
the  problem,  it  is  not  possible  to  make  a  definitive  comparison  of  the 
relative  efficiency  of  the  subject  codes. 

In  the  next  section,  we  analyze  the  reduced  model  for  sensitivity 
information  and  optimal  value  parametric  bounds. 


5.  Optimal  Value  Function  Sensitivities  and  Bounds 


5.1  Optimal  value  function  sensitivities 

To  conduct  sensitivity  and  bound  calculations  with  respect  to 
the  various  parameters  involved  in  the  model,  the  coefficient  c^  (1*1,  20) 

must  be  expressed  as  a  function  of  the  model  parameters,  rather  than  being 
considered  independent  and  assigned  fixed  values  (as  was  done  in  defining 
the  data  of  the  problem  solved  in  Section  4) .  To  accomplish  this,  the 
formulas  taken  from  [4]  and  given  in  Section  3  for  c^  (i*l,  20)  were 

coded.  Using  the  values  for  the  various  problem  parameters,  given  in 
1 3 1  and  [15]  and  shown  in  Table  1,  the  following  values  were  calculated  for 
the  components  of  the  parameter  -  dependent  coefficient  vector  c*('c^, . . .  .c^q)  . 

2.477010  *  10* 

1.028863  *  10 

8.166663  *  10 

1.0 

.08333296 


C  L  = 

1270.5 

C6 

=  .2397342 

C11  = 

6.213914 

*  10~8 

C.  m 

16 

C2  = 

.01229334 

C7 

=  1.0 

C12  = 

50.0 

C17  " 

-9 

C3  = 

2.419023 

C8 

=  .6200489 

C13  " 

4.739336 

*  10 

C18  * 

-3 

-6 

c4  = 

77169.87 

C9 

=  1.107162  *  10 

C14  = 

7.312399 

*  10 

f*  ■ 

C19 

C5  * 

888.7634 

C10 

-  2.18723  *  1018 

C15  “ 

1.2577  * 

1C'5 

c  K 

20 

-3 


-3 
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Table  1 
Problem  Data 
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except  i  or  c  ,  the  components  of  c  are  In  close  agreement  with  those 

used  in  [  1]  and  given  in  Section  4.  (We  note  that  in  113],  the  value  ol 
c  conforms  exactly  to  the  above  value,  which  is  i/2  of  the  value  reported 

in  i 3] .  This  suggests  chat  the  vaLue  of  c  derived  here  is  correct.)  The 

results  reported  in  the  remainder  of  the  paper  are  consequently  based  on 
Lae so  components  of  c  . 

Table  2  depicts  the  sensitivity  of  the  optimal  total  cost  with 
respect  to  numerous  parameters  of  the  model.  As  shown  in  this  table, 
we  do  not  conduct  sensitivity  analysis  with  respect  to  all  of  the  model 
data  (parameters)  listed  in  Table  1.  This  is  because  some  data  cannot  be 
meaningfully  altered.  For  example,  the  gravitational  constant  g  is 

assumed  known  and  fixed.  The  specific  heat  c  ,  Che  viscosity  p  ,  etc., 

P 

are  internal  properties  of  the  fluid  and  cannot  freely  be  altered. 

As  elaborated  in  our  earlier  papers,  [8],  [9],  the  sensitivity 
measure  df*(i  )/l)ci  corresponds  to  the  rate  of  change  of  the  optimal 

value  function  f*(t.)  with  respect  to  t.  ,  at  a  given  value  of  £  . 

It  would  seem  that  the  larger  the  absolute  value  of  df*(t;)/3e  , 

the  greater  the  sensitivity  of  f*(i.)  with  respect  to  e.  .  However, 

it  would  be  erroneous  to  conclude  this  and  misleading  to  use  the  informa¬ 
tion  without  further  analysis,  as  a  guideline  for  making  and  implementing 
appropriate  decisions.  This  is,  in  part,  because  this  rate  of  change  is 
Local  and  does  not  necessarily  reflect  the  change  that  can  result  from 
finite  data  changes.  Hor  does  this  measure  reflect  the  likelihood  of  a 
given  finite  change.  For  such  reasons,  it  is  imperative  to  calculate 
some  measures  of  stability  in  the  large,  to  estimate  or  at  least  contain 
the  ctfects  of  finite  data  changes  on  the  optimal  solution.  Moreover, 
sensitivities  are  often  interpreted  as  the  "change  ir.  f*(t;)  per  unit 

change  in  t  ^  "  .  Such  interpretations  may  also  prove  to  be  misleading 

because  they  tacitly  suggest  that  a  unit  change  is  meaningful  and  that  it 
might  be  realized  in  a  given  parameter  .  However,  a  one  unit  change 

lor  one  parameter  may  be  too  large  to  be  practically  feasible,  while  for 
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Optimal  Value  Function  Sensitivity  Results 

—  —  (5) 

Parameter  Scaled 

Value  Sensitivity  Sensitivity 

e±  3f*(e)/9e.  |  .01e±)  Of*(e)/3ei) 


-159.9 

22.577.27 

36,101.05 

•  32.08 

-31,283.56 

-10,035.77 

1.1087 

-.3294094  *  108 

-365,216.2 

.7127 

341,840.9 

2,436.30 

92340 

8.288348 

7,653.46 

27.444 

22,112.94 

6,068.67 

1.2 

33,954.16 

407.45 

1.7599 

1,770.420 

31,157.62 

1.0 

-2,453,125. 

-24,531.25 

1.1407  *  10-3 

-.2583156  *  1010 

-29,466.06 

1.2716  *  10-3 

.2680632  *  108 

341.14 

3.0122 

3,316,445. 

99,897.95 

3.4851  *  10-5 

.1728371  *  1011 

6,023.55 

1.5567  *  104 

28.58427 

4,449.71 

1. 

4,710,256. 

47,102.56 

8250. 

80.40698 

7,376.07 

.00067 

.4588196  *  1010 

30,740.91 

5.0 

55,620.45 

2,781.02 

.06408 

3,797,937. 

154.96 

continued 
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Table  2  (continued) 

Optimal  Value  Function  Sensitivity  Results 


(1) 

Parameter  No. 

1 

(2) 

Parameter 

'i 

(3) 

Parameter 

Value 

Ei 

(4) 

Sensitivity 

9f*(c)/Df.1 

(5) 

Scaled 

Sensitivity 

|.01c1|(3f*(c)/3e1) 

20 

P 

c 

.154 

6,609,746. 

10,179.01 

21 

T 

c 

50.0 

17,350.95 

8,675.47 

22 

1 1 

6337.86 

4,850.329 

307,407.06 

23 

w 

1.055  *  108 

-.0175986 

18,566.52 

24 

+  r>G 

.988 

.3107304  *  108 

307,001.64 

25 

Pmi 

1.25(»  — 

32,596.29 

407.453 

26 

Qa,R 

2.713  *  106 

12.97611 

352,041.86 

27 

h2,R 

989.1 

31,283.35 

309,423.61 

28 

fM,R 

.7217 

-684,874.7 

-  4,942.74 

t  Candidate  for  parametric  bound  calculation  on  f*(e)  . 
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another  it  may  be  too  small  to  be  numerically  significant.  Parameters 
n  .  the  turbine  efficiency,  and  W  ,  the  water  (low  rate,  in  Table  2 

(j 

provide  good  exampLes  of  boLh  Instances. 

With  these  points  in  mind,  and  assuming  that  the  sensitivities 
are  acceptable  in  a  small  (finite)  neighborhood  of  the  problem  data, 
we  prefer  to  base  our  judgments  on  "scaled"  sensitivities  such  as  those 
given  in  Column  (5)  of  Table  2,  which  correspond  to  the  estimated  changes 
resulting  from  a  one  percent  change  in  a  given  parameter  value.  Column 
6  of  this  table  ranks  the  above  scaled  sensitivity  measures,  in  decreasing 
order  of  their  absolute  values.  The  importance  of  providing  the  kind  of 
information  depicted  in  Table  2  for  real  world  systems  seems  obvious  and 
cannot  be  overstated.  For  example,  it  provides  guidelines  for  setting 
priorities  on  the  relative  importance  of  the  parameters,  and  suggests 
actions  that  can  be  taken  to  obtain  optimum  marginal  improvements  in 
the  optimum  performance. 


5.2  Optimal  value  function  bounds 

here,  we  apply  the  bound  calculation  technique  discussed  in 
Section  2,  to  derive  piecewise  linear  parametric  upper  and  lower  bounds 
on  the  optimal  total  cost  of  the  steam  turbine  power  plant  under  study. 


The  parameters  marked  by  (t)  in  Table  2  are  retabulated  in 
table  1  and  are  selected  for  bound  calculation.  Notice  that  parameter 
a.j  is,  in  fact,  a  right-hand  side  parameter  of  constraint  g^  .  Recall 

that 

*5  =  C16  X8_1  +  C17  X2  -  1 


where 


'16 


a,R 


and  c.,  = 

1  i 


Thus  g,.  can  be  written  as 


Qa,R  x8 


-1 


-  b3  X2  i  *3 


Thus,  optimal  value  function  of  the  model  under  study  is  a  convex  non¬ 
increasing  function  of  a^  . 
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Table  3 

Parameters  Selected  For  Bound  Calculation 


Pa rameter 

Parameter 

Parameter 

Perturbation 

No. 

Name 

Value 

3 

°3 

1.1087 

.2527 

16 

c 

8250. 

-2062.5 

an 

18 

ct 

5.0 

-1.25 

22 

t 

6337.86 

-1584.465 

24 

nc 

.988 

-.247 

25 

p  ,(*  — ) 
mi  np 

1.25  (-  ~) 

-.125 
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Also,  parameters  e  ,  r  ,  t  and  n  .  taken  one  at  a  time,  appear  In 

the  objective  function  cue! f  If  lent a  linearly  and  do  not  appear  In  the 
constraints.  Thus,  the  optimal  annual  cost  of  the  plant  under  considera¬ 
tion  is  a  concave  nondecreasing  function  of  each  of  these  parameters. 

This  means  that  our  bound  calculation  procedure  can  readily  be  applied  to 
these  parameters.  Notice  that  parameter  np  (pump  efficiency)  does  not 

appear  linearly  in  the  objective  function  coefficients.  For  this  reason, 
we  replace  — - -  by  P  .  .  This  enables  us  to  calculate  optimal  value 

Hp  mi 

bounds  relative  to  changes  in  P  .  and,  hence,  relative  to  changes  in  n  • 

mi  r 

The  last  column  in  Table  3  shows  the  magnitude  of  the  perturbations 

allowed  in  the  respective  parameter  value,  each  taken  in  the  direction 

that  decreases  the  optimal  value  function,  as  determined  by  the  sensitivity 

information  In  Table  2.  In  the  following,  we  list  the  resulting  parametric 

bounds  on  the  optimal  value  over  the  given  range  of  the  parameter,  obtained 

by  our  latest  version  of  the  computer  programming  code  SENSUMT  [10]. 


•  Hounds  w.r.l.  a,  ,  1.1087  >-  a..  •  j.3614 

>  a 

Upper  bound  V  =  -  .1792469  *  108  «3  +  .5163132  *  10* 

Lower  bound  K  =  fiax  {-  .3294094  *  10°  a3  +  .6827984  *  108, 

-  .6760806  *  107  a  4  . 3b43/82  *  1C8} 

Figure  2  depicts  the  resulting  computer  output  lor  this  and  supplemental 
bound  calculation  for  a3  . 

•  liounds  v.r.t.  c  ,  6187.5  <  c  <  S25G  . 

an  an  - 
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Optimal  Value  Function  Bounds  When  Par(  3)  is  Perturbed 


Point  1  (Unperturbed  Solution)  : 

Par(  3)  «  0. 1108700D  01  F(Par(  3))  =  0.3175823D  08 

Point  2  (Perturbed  Solution)  : 

Par (  3))  =  0. 1361400D  01  F(Par(  3))  =  0.2722866D  08 

Line  passing  through  Points  1  and  2  and  overestimating  f* 

F  =  -0. 17924691)  08  *  Par(  3)  +  0.5163J32D  08 
Line  underestimating  f*  at  Point  1 

F  =  -0.3294094D  08  *  Par(  3)  +  0.6827984D  08 

Line  underestimating  f*  at  Point  2 


F  =  -0. 6760806D  07  *  Par(  3)  +  0.3643282D  08 
Quadratic  Estimation  of  f*  through  Points  1  and  2 

F  =  0.441784  ID  08  *  Par(  3)**2  -  0.1270497D  09  *  Par(  3)  +  0.1183135D  09 
1*  bound  evaluation  at  ten  equidistant  points  between  Points  1  and  2 


Par(  3)  Lower  Bound  Upper  Bound  Quad  Estimate 


1.109 

0. 3175823D 

08 

0. 3175823D 

08 

0. 3175823D 

08 

1.134 

0. 3092581D 

08 

0. 3130527D 

08 

0. 3105137D 

08 

1.159 

0. 3009339D 

08 

0.3085231D 

08 

0. 3040093D 

08 

1.185 

0.2926097D 

08 

0.3039935D 

08 

0. 2980692D 

08 

1.210 

0. 2842856D 

08 

0. 2994640D 

08 

0.2926933D 

08 

1.235 

0.2808288D 

08 

0. 2949344D 

08 

0.2878816D 

08 

1.260 

0. 2791204D 

08 

0.2904048D 

08 

0.2836342D 

08 

1.286 

0. 2774119D 

08 

0. 2858753D 

08 

0.2799509D 

08 

1.311 

0.2757035D 

08 

0. 2813457D 

08 

0. 2768319D 

08 

1.336 

0.2739950D 

08 

0. 2768161D 

08 

0. 2742771D 

08 

1.361 

0.2722866D 

08 

0. 2722866D 

08 

0. 2722866D 

08 

Figure  2.  Parametric  bounds  on  f*(a^) 
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Upper  bound 

Lower  bound 

.  Bounds  w.r.t 

Upper  bound 

Lower  bound 

Figure  3  depicts  the 
respect  to  c  . 

.  Bounds  w.r.t 
Upper  bound 

Lower  bound 

.  Bounds  w.r.t 

Upper  bound 

Lower  bound 


F  -  min  {89.40698  c  +  .3102062  *  10®  , 

un 

100.5829  c  +  .3094000  *  108} 
an 

F  =  94.96585  c  +  .3097476  *  10® 

an 

.  ct»  3.75  <  <  5.0 

F  *  min  {55620.45  c  +  .3148012  *  108  , 

60878.7  ct  +  .3145717  *  108} 

F  =  58209.10  ct  +  .3146718  *  108 

resulting  computer  output  for  bound  calculation  with 

t,  4753.395  <  t  <  6337.86 

F  =  min  {4850.329  t  +  .1017524  *  107  , 

4860.131  t  +  .8543850  *  106} 

F  *  4863.883  t  +  .9316157  *  106 

nG  ,  .741  <  nG  5  -988 

F  =  min  {.3107304  *  108  +  .1058064  *  107  , 

G 

.3126812  *  108  nG  +  .8912627  *  107} 

F  -  .3116310  *  10®  nG  +  .9690873  *  106 
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Optimal  Value  Function  Bounds  WUen  Far (18)  is  Perturbed 


Point  1  (Unperturbed  Solution)  : 

Par(  18)  -  0.50000001)  01  F(Par(18))  -  0.3175823D  08 

Point  2  (Perturbed  Solution)  : 

Par(18))  -  0.37 500 00D  01  F(Par(18))  -  0.3168546D  08 

Line  passing  through  Points  1  and  2  and  underestimating  f* 

F  -  0. 5820910D  05  *  Par (18)  +  0.3146718D  08 
Line  overestimating  f*  at  Point  1 

F  -  0. 5562045D  05  *  Par(18)  +  0.3148012D  08 
Line  overestimating  f*  at  Point  2 

F  -  0.6087870D  05  *  Par(18)  +  0.3145717D  08 
Quadratic  estimation  of  f*  through  Points  1  and  2 

F  -  -0.2070913D  04  *  Par(18)**2  +  0.7632959D  05  *  Par(18)  +  0.3142835D  08 
f*  bound  evaluation  at  ten  equidistant  points  between  Points  1  and  2 


Par (18) 

Lower  Bound 

Upper  Bound 

Quad  Estimate 

5.000 

0. 3175823D  08 

0.3175823D  08 

0.3175823D  08 

4.875 

0. 3175095D  08 

0.3175127D  08 

0.3175124D  08 

4.750 

0. 3174367D  08 

0.3174432D  08 

0. 3174419D  08 

4.625 

0. 3173640D  08 

0.3173737D  08 

0.3173708D  08 

4.500 

0.3172912D  08 

0. 317304 2D  08 

0. 3172990D  08 

4.375 

0. 3172184D  08 

0. 3172346D  08 

0. 3172265D  08 

4.250 

0. 3171457D  08 

0.3171590D  08 

0. 3171535D  08 

4.125 

0.3170729D  08 

0.3170829D  06 

0. 3170797D  08 

4.000 

0.3170002D  08 

0.3170068D  08 

0. 3170Q53D  08 

3.875 

0. 3169274D  08 

0.3169307D  08 

0.3169303D  08 

3.750 

0.3168546D  08 

0.31 68 54 6D  08 

0.3168546D  08 

Figure  3.  Parametric  bounds  on  f*(ct> 
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Bounds  on  P  ,  = 

1 

-  ,  1.125  < 

1 

< 

mi 

T. 

i 

nP 

Upper  bound  F  = 

min 

{.3259629  * 

io5 

1 

n 

.3536232  * 

105 

l 

n 


Lower  bound  F  =  .3355953  *  105  — - —  + 

nP 

These  bounds  are  piecewise  Linear  as  a  function  of 
nonlinear  as  a  function  of  np  • 


1.25  (.8  <  np  <  .889) 

—  +  .3171748  *  108  , 

P 

—  +  .3171425  *  103} 

P 

.3171628  *  108 


P  .  ,  hence  they  are 
mi 
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